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Abstract 



Feedback control laws have been traditionally implemented in a periodic fashion on digital hardware. 
Although periodicity simplifies the analysis of the mismatch between the control design and its digital 
implementation, it also leads to conservative usage of resources such as CPU utilization in the case of 
embedded control. We present a novel technique that abandons the periodicity assumption by using 
the current state of the plant to decide the next time instant in which the state should be measured, 
the control law computed, and the actuators updated. This technique, termed self-triggered control, is 
developed for two classes of nonlinear control systems, namely, state-dependent homogeneous systems 
and polynomial systems. The wide applicability of the proposed results is illustrated in two well known 
physical examples: a jet engine compressor and the rigid body. 



1 Introduction 



Feedback control laws are mainly implemented nowadays in digital platforms since microprocessors offer 
many advantages with respect to analog platforms. Despite its advantages, digital implementations raise 
several design difficulties; one of the most relevant problems consists in determining how frequently the 
controller needs to be executed so that a desired performance is achieved. Traditionally, the controller is 
executed periodically every T units of time, since this facilitates the analysis of the mismatch between the 
customary continuous design and the digital implementation. However, it seems unnatural to update the 
signals of interest in a periodic fashion, especially for nonlinear control systems: control tasks should be 
executed only when something significant happens in the process to be controlled. Moreover, the choice 
of the period T is based on a worst-case scenario (to guarantee performance for all possible operation 
points), and hence the control task is executed at the same rate regardlessly of the state of the plant. 
Therefore, the periodicity assumption leads to inefficient implementations in terms of processor usage, 
communication bandwidth, energy... For instance, in the context of embedded systems, the CPU time 
is shared between different tasks. Hence, executing the control task when the states barely change is a 
waste of computational resources. In the context of distributed control systems, not only the processor 
time is a scarce resource but also the available communication bandwidth. Thus communication should 
only occur when relevant information needs to be transmitted from the sensor to the controller and/or 
from the controller to the actuator. 

There is a vast literature about periodic implementations of control laws; still, there are many open 
points, especially for nonlinear systems. We refer to |HL05] for a good introduction to the subject for 
linear systems. Due to the lack of a deep understanding of the effects of digital implementations, ad 
hoc rules are commonly used to determine stabilizing periods (for instance, 20 times the time constant 
of the dominant pole |AW90j ) . Most of the results relating stability and sampling periods for linear 
systems are based on the construction of an equivalent discrete-time model. Nonlinear systems, in 
general, cannot be discretized exactly in closed form. Hence, a common approach consists in finding 
approximate discrete-time models |NTK99] . and then carry out the study and design for this set of 
equations [LNA05] . 

To overcome the drawbacks of the periodic paradigm, several researchers ( |Arz99 , |AB02j . |San06] . 



[Tab07]) suggested the idea of event-triggered control. Under this paradigm, the controller execution 
is triggered according to the state of the plant. The event-triggered technique reduces resource usage 
and provides a high degree of robustness (since the system is measured continuously). Unfortunately, 
in many cases it requires dedicated hardware to monitor the plant permanently, not available in most 
general purpose devices. 

In this paper we propose to take advantage of the event-triggered technique without resorting to extra 
hardware. In most feedback laws, the state of the plant has to be measured (or estimated) to compute 
the next value of the controller; hence, this information could also be used to decide when the controller 
needs to be executed again. This technique is known as self-triggered control, since the controller 
decides its next execution time. This technique could be regarded as a way to introduce feedback in 
the triggering process, in comparison to the open loop structure in periodic implementations. A first 
attempt to explore self-trigger models for linear systems was developed in [VFM03 , by discretizing 
the plant, and in (WL08 for linear Tioo controllers. No studies have been carried so far for nonlinear 
systems, where aperiodicity is expected to be even more beneficial since the behaviour of a nonlinear 
dynamical system highly depends on the operation point. 

The contributions of this paper are twofold: 

1. Under self-trigger implementations, the execution times for a control law are defined by a simple 
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self-trigger condition (since it has to be computed online) that depends on the dynamics of the 
system, the desired performance, and the current measurement of the state. We develop self- 
trigger conditions for two classes of systems, namely, state-dependent homogeneous systems and 
polynomial systems. 

2. The developed self-trigger conditions are but a facet of the inherent tradeoff between the amount 
of resources alloted to a controller and the achieved performance. On our route to develop these 
conditions we develop the necessary theory required to further the understanding of the real-time 
requirements of control systems. 

The results in this paper can thus be seen as a contribution to the broad area of control under 
information constraints. Many authors have recently studied other aspects of this problem, such as 
space quantization [BLOOj, delays in network control systems [LMT03 or minimum attention control 
design |Bro97] , to cite a few. 

The rest of the paper is organized as follows. We first review in section III an event-trigger condition 
that guarantees stability under sample-and-hold implementations, previously studied in |Tab07 . In 



section IV we define homogeneous control systems and we analyse the properties of their trajectories. 
Such properties can be exploited to derive a simple self-trigger condition based on the norm of the 
current state. This idea is generalized by introducing the notion of state-dependent homogeneity, that 
allows us to establish a more general triggering condition based on the full state of the plant. Finally, 
in section V the previous results are utilized to develop self-trigger conditions for (^-related vector fields 
and, in particular, for polynomial control systems. To illustrate the wide applicability of the results we 
select two very different examples from the literature: the control of a jet engine compressor and the 
control of a rigid body. The appendix includes the proofs for all the results appearing in the paper. 



2 Notation and input-to-state stability 
2.1 Notation 

We shall use the notation \x\ to denote the Euclidean norm of an element x £ W 1 . A continuous function 
a : [0, a[— > Wq, a > 0, is said to be of class /C if it is strictly increasing and a(0) = 0. It is said to be of 
class /Cqo if o, = oo and a{r) — * oo as r — > oo. 

A function is said to be of class C°° or smooth if it can be differentiated infinitely many times. All 
the objects in this paper are considered to be smooth unless otherwise stated. 

Given vector fields X and Y in an n-dimensional manifold M, we let [X,Y] denote their Lie product 
which, in local coordinates x = (x%,X2, ...,x n ), we take as j^(x)Y(x) — ^(x)X(x). Let <p : M — ► N be 
a map, and let X and Y be vector fields in manifolds M and N, respectively. The differential of a map 
ip is denoted by Tip, and whenever we need to emphasize the point x where the differential is evaluated 
we write T x ip. We call X and Y 99-related if the following holds: 

Tip X = Y o <p 

If <p has a smooth inverse, the pullback of the vector field Y under the mapping tp can be defined as: 

ip*Y = (Tip)- x Yoip 

The notations ip(t, •) and ipt(') will be used interchangeably to denote a map i/):ixM-> M. Finally, 
we use e s to represent the exponential of s G R. 
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2.2 Input-to-state stability 

We consider a control system: 

x = /(x,u), 



x(t) g M n , u(t) g M m ,te 



(2.1) 



We use x G M n to denote the state of the control system, x(i) to denote a solution of (2.1 ), and u(t) for 
the trajectory of the input. Whenever it is needed to emphasize the initial condition for the trajectory 
we write x(t,xo), where x(0,xq) = xq. We will not resort to the standard definition of input-to-state 
stable (ISS) (see |Son05 for details) in this paper but rather to the following characterization. 

Definition 2.1 A smooth function V : 1" — ► IR^~ is said to be an ISS Lyapunov function for the 
system (2.1) if there exist class /Coo functions a, a, a and 7 satisfying: 



a(\x\) < V(x) < a{\x\) 
fK/(x,u)<-a(M) + 7 (M) 



(2.2) 
(2.3) 



The system (2.1) is said to be ISS with respect to the input u if and only if there exists an ISS lyapunov 
function for (2.1). 



3 Event-triggered stabilization of nonlinear control systems 

Even though most of the current controllers are implemented on digital platforms, the techniques used to 
design these controllers ignore the implementation effects, such as delays or space and time quantization. 
In this paper we focus on the time quantization aspects. We start analyzing the behaviour of the plant 
under the event-triggered implementation introduced in |Tab07] and reviewed in this section. Consider 



the control system described in (2.1) for which a feedback controller u = k(x) has been designed. The 
implementation of such feedback law on an embedded processor is typically done by sampling the state 
at time instants t{, computing u(ij) = A;(x(tj)) and updating the actuator values at time instants tj + Aj. 
For ease of exposition, it is assumed that Aj = for all z G N, that is, every time the state is sampled the 
controller is computed without dela>Q Furthermore, the sequence of times {tijieN is typically periodic 
meaning that ti + \ — ti = T for any i 6 N, where T > is the period. In this paper the periodicity 
assumption is dropped in favour of self-trigger implementations. We would like to identify a class of 
aperiodic executions that guarantees stability and desired performance while reducing the number of 
executions. In other words, we would like to determine when it is indeed necessary to execute the 
control law in order to achieve the desired performance. To derive a stabilizing triggering rule, we look 
at the mismatch between the current state and the sampled state. Similar results could be developed by 
looking at the mismatch of the input of the plant under analog and digital implementations. Since both 
analysis are analogous from a conceptual point of view, in this paper we just focus on the state-based 
approach for ease of exposition. 

We define the measurement error e as the difference between the last measured state and the current 
value of the state: 

e(t) = xfe) - x(i) for t G [t h t i+1 [ (3.1) 
With this definition, the closed loop x = /(x, fc(x(tj))) becomes: 

i=/(x,fc(x + e)) (3.2) 
lr The results can be generalized for nonzero A by following the procedure described in [Tab07]. 
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Let the control law u = k{x) render the system ISS with respect to the measurement error e. Under this 
assumption, there exists a Lyapunov function V for the system that satisfies the following inequality: 



V<-a{\x\)+ 1 {\e\) 



(3.3) 



where a and 7 are /Coo functions. Stability of the closed loop system (3.2) can be guaranteed if the 
error is restricted to satisfy: 

l(\e\)<aa(\x\), a > (3.4) 
since the dynamics of V will be bounded by: 

V <{(X-l)a{\x\) (3.5) 

guaranteeing that V decreases provided that a < 1. If a~ l and 7 are Lipschitz continuous on compacts, 



inequality (3.4) is implied by the simpler inequality: 



b\e\ < aa\x\ 



(3.6) 



for a and b appropriately chosen according to the Lipschitz constants of a 1 and 7. Inequality (3.6) 
can be enforced by executing the control task whenever: 



a 

cr-\x\ 
b 



(3.7) 



Upon the execution of the control task, the state is measured and the error becomes 0, since x(ij) = x(i) 
implies e(t) = x(tj) — x(t) = 0. An event-triggered implementation based on this equality would require 



testing (3.7) frequently. Unless this testing process is implemented in hardware, one might run the 



risk of consuming the processor time testing (3.7). To overcome this drawback, we opt for self-trigger 
strategies, where the current state measurement x(ij) is used to determine the next execution time ti + \ 
for the control task. 

to evolve from 2 to 



The inter-execution time implicitly defined by (3.7) is the time it takes for M 



u 



t. Since this time depends on the last sample of the system x(ij), it will be denoted as r(j;(tjjj. In 
order to derive self-trigger conditions based on the sampled state, we study in the next sections the 
evolution of r as a function of x{ti). This evolution can be described by analyzing the trajectories for 
different types of nonlinear systems. For simplicity of presentation, we start considering homogeneous 
systems of constant degree. 



4 Homogeneous control systems 

Homogeneous vector fields are vector fields possessing a symmetry with respect to a family of dilations. 
They appear as local approximations for general nonlinear systems |Hcr91J since an analytic function 
can always be decomposed in an infinite sum of homogeneous functions. Moreover, many physical 
systems can be described as homogeneous systems (see |Bai80| for some examples). 

4.1 Definitions 

To define homogeneity we first review the notion of dilation. 

2 Recall that at the execution instant t — ti we have e(t) — x(ti) — x(t) = and thus W- = 0. 
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Definition 4.1 Given an n-tuple r=(n, ..,r n ) £ (IR^~) n , i/ie dilation map 5 r x : W 1 — > M n is defined by: 

5 r x (x) = (X^x 1 ,...,X r "x n ), A>0 (4.1) 
^4 dilation is called standard if r{ = /or all i,j £ {1, . . . ,n}. 

Definition 4.2 >1 function f : M n — > M n is called homogeneous of degree £ if for all A > 0, t/iere exist 
r = (ri,..,r n ) £ swc/i t/iat; 

/(<^(x)) = A^(/(x)) (4.2) 

where C > — minj r i . 

With this definition, we see that linear functions are homogeneous of degree C = with respect to the 
standard dilation. Consider now a differential equation: 

x = /(x) (4.3) 

whose right-hand side is homogeneous of degree Q. Then, the solution 4>(t, 0o) satisfies: 

0(i,^(xo)) = <W(A^,z o ) ( 44 ) 



For a proof see |Kaw95] or the proof of Theorem 4.5 in this paper that generalizes this fact. Equality (|4.4 



lies at the heart of the results presented in this paper. 

4.2 Scaling laws for the inter-execution time of homogeneous systems 

In this section we consider homogeneous control systems with respect to the standard dilation, since 
the general case can be reduced to this one, as explained in [GriiOO . Using the commutative property 



of homogeneous systems expressed by equality (4.4), we can derive the following scaling law for the 



inter-execution times under the event-triggered policy described in Section [3] 

Theorem 4.3 Let x = /(x, u) be a control system for which a feedback control law u = k(x), rendering 
the closed loop homogeneous of degree £ with respect to the standard dilation, has been designed. The 
inter- execution times r : W 1 — > R implicitly defined by the execution rule |e| = c|x| with c > scale 
according to: 

t(S t x (x)) = A" c r(x), A > (4.5) 
where x € M n represents any point in the state space. 

Similar results were obtained in |Tun05 but with a different objective (finite-time stability and finite 
escape time for homogeneous systems). This theorem relates the inter-execution time at a point x with 
the inter-execution time at any other point lying along the same homogeneous ray. That is, once the 
time is known for just one initial condition x, we can infer the times for all initial conditions of the type 
Ax, for any A > 0. In the next section these ideas are extended by allowing the degree of homogeneity 
C to be a function of the state. 
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4.3 State-dependent homogeneity and space-time dilations 



In order to generalize the classical notion of homogeneity presented in section [4] we resort to the 
coordinate-free geometric notion of homogeneity introduced in |Kaw95] . 

Definition 4.4 Let D : M — > TM be a (dilation) vector field such that x = —D(x) is globally asymp- 
totically stable. A vector field X : M — > TM is called homogeneous of degree ( with respect to the vector 
field D if it satisfies the following relation: 

[D,X] = CX (4.6) 
The trajectories tp : K x M — > M of D are called the homogeneous rays of the system. 
To recover the classical notion of homogeneity, we set: 

n d 

J D = Vr i x i — , ri >0 (4.7) 

r-f oxi 

i=i 

where { ^r} i _ 1 n is the canonical basis for the tangent bundle TM. By substituting this expression for 
D in ( |4.6| , the classical notion of homogeneity (equation ( |4.2[ ) ) can be recovered as explained in |Arn92| . 
Under this framework, we can see that the definition of section [4] just corresponds to a particular choice 



of the vector field D. In this geometric context, (4.4) can be generalized to the following theorem. 



Theorem 4.5 Let X and D be vector fields on a manifold M, giving rise to flows (ft : R x M — > M and 
ifj : M x M — > M respectively. The vector field X is homogeneous of degree ( with respect to D if and 
only if: 

cftt o i> s = tj} s o (ft e(st , a, t € R (4.8) 

This theorem implies that the flows of homogeneous vector fields commute in a particular way: applying 
the flow ip before the flow (ft to a point x £ M is equivalent to applying the flow (ft before ip but with 



a scaling in time, given by e^ s , where £ is the degree of homogeneity. For D as defined in (4.7), 
equation (4.8) simplifies to (4.4), since the flow ift s becomes the dilation map 5X, for A = e s . In this 



paper we generalize even further Definition 4.4 to the following state-dependent notion of homogeneity. 



Definition 4.6 Let D : M — > TM be a (dilation) vector field such that x = —D(x) is globally asymp- 
totically stable. A vector field X : M — > TM is called homogeneous with degree function £ : M — * M 
with respect to the vector field D : M — » TM if it satisfies the following relation: 

[D,X] = £X (4.9) 

Example 4.7 Let the vector fields D and X be: 

d d 

D = ax\— h OJX27: — j a>0 

oxi 0x2 

X = {- Xl - x \)JL + {- X2 - x \ X2 )A- (4.10) 
ox\ 0x2 



The vector fields D and X do not satisfy (4-6) for any value of (, but they satisfy |^.g[ ) for the following 
function 

3xj + 1 
x\ + 1 
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As for the standard homogeneity case, a relation between trajectories can also be derived. 

Theorem 4.8 Let X and D be vector fields on a manifold M, giving rise to flows <p :lx M — > M and 
^:RxM->M respectively. The vector field X is homogeneous with degree function £ : M — > R witt 
respect to D if and only if: 

4>t ° ips = i>s ° 4> e p(s) t s,t£R (4.11) 
p( s ) = ^o^dr (4.12) 



./o 

The vector field D generates a flow ip that can be considered as a spatial dilation operator acting 
on points in M. This operator determines how times are scaled in the flow of 4>. That is, given a point 
x in a manifold M, applying this spatial operator tf> to x entails a scaling in time e p ^ in the flow of a 
vector field X. 



Remark 4.9 In the case of the standard homogeneity, we have £(x) = ( and thus equation (4-11) 



becomes (4-8) since p is linear in time, p(s) = (s. Furthermore, for linear systems, the degree of 



homogeneity £ is 0, therefore p(s) = and the vector fields X and D commute, that is: 

4> t o^ s = ip s o (p t (4.13) 



Example 4.10 To illustrate the previous theorem, we recover Example 4-7 The flow for the dilation 
vector field D is: 

il> s (x)=e as x (4.14) 



where x = (xi, X2) T ■ Hence the flow of the vector field X satisfies equation (4-H )'■ 

d> t (e as x) = e a ^ eP{s)t (x) (4.15) 

with: 

We see that the dilation in space e as induces a scaling in time e p ^ in the trajectory <f>. 

Using these theorems, a self-trigger condition for state-dependent homogeneous systems can developed. 

4.4 Scaling laws for inter-execution times of state-dependent homogeneous systems 

For simplicity, we use the standard dilation vector field: 

n d 

D = Trxi— (4.17) 

r-f oxi 

i=i 

The flow of this vector field (i.e., the homogeneous ray) is: 

ipi(s,x) = e rs Xi i = l,...,n x = (xi, . . . , x n ) T (4.18) 



Using the commutative property (4.11) of homogeneous flows, we can derive a generalized scaling law 
for the inter-execution times. 
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Theorem 4.11 Let x = /(x, u) be a control system for which a feedback control law u = k{x) has been 
designed, rendering the closed loop homogeneous with degree function £ with respect to the standard 
dilation vector field. The inter- execution times r : W 1 — > R implicitly defined by the execution rule 
\e\ = c\x\, with c > 0, scale according to: 

r(e rs x) = e~ p ^r(x) Vs G R, r > 

with: 



(4.19) 



£(e rv x)dv 



(4.20) 



o 



and where x G M. n represents any point in the state space. 

Remark 4.12 For linear systems, the degree function £ is 0, therefore p(s) = and e^ p ^ = 1. That 
is, for linear systems the inter- execution times remain constant as we move along homogeneous rays. 



Remark 4.13 When the degree function £ is a constant (Definition 4-4h ^ e scaling in time depends 
entirely on the value of s. That is, to see how times scale it is enough to determine the position of 
the state x on a homogeneous ray, but not on which particular ray x lies, since time in all rays scale 
in the same way. In other words, the scaling in time is just a function of the norm of the state, and 
therefore (4-19) simplifies to (4-5), for A = e s and p(s) = £s. On the other hand, for state-dependent 
homogeneity the scaling in time is determined by the ray where x lies and the position of x in that ray, 
that is, p is a function of s and x, since time in each homogeneous ray scales in a different manner. 

Theorem 14.111 allows us to use the estimate of the inter-execution times at some x in order to 
determine the inter-execution times for the whole ray through x. Therefore, it is enough to find estimates 
of these times on any n — 1 sphere, and then extend the estimates along homogeneous rays. Moreover, 
since a lower bound r* for the inter-execution times can be easily computed for linear system^] a n — 1 
sphere can always be chosen where a linear over-approximation for the control system can be obtained. 
To do so, we rewrite equation (4.9) in local coordinates. A homogeneous function g of degree £(x) — 1 
satisfies: 

n ft 

+ n) g {x) = riXi ]^r. ( 4 - 21 

i=l 



Hence, for the closed loop system x 
\x\ and lei: 



/(x, e) = /(x, k(x + e)) we can find a bound for |/(x, e)| linear in 



\f(x, e)| = \H(x, e)x + G(x, e)e| 

< \H(x,e)\\x\ + |G(x,e)||e| 
<\H(x* a ,e* a )\\x\ + \G(x* b ,el)\\e\ 



where: 



H(x,e) 



G{x,e) 



ri d/i 

ri dfi 
£+ri dei 

r n dfn 
- £+r„ dei 



ri dfi 
£+ri dx n 

r n df n 
£+r„ dx n 

ri dfi 
f+ri de n 

r n dfn 
i+r n de n 



(4.22) 



s Scc TabOr for one possible method to find such r* 
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and (x*,e*) and (x* b ,e* b ) are such that: 

\H(x,e)\<\H(x* a ,e* a )\, \G(x, e)\ < \G(x* b , e* b )\ 

for all (x, e) in a neighbourhood 0, around the origin. So given this set O we can find where the norm 
of these matrices H and G attain its maximum values and then work with the following linear model 
for the sole purpose of computing r* : 

± = H(x* a ,e* a )x+G(x* b ,e* b )e (4.23) 

It is important to emphasize that we are not trying to find a linearized model, as it would not guarantee 
stability for the original nonlinear system. To summarize, the computation of the self-triggered execution 
strategy is made in 4 steps: 

1. Define an invariant set O x C M. n around the equilibrium point, for instance a level set of the 



Lyapunov function. The execution rule (3.7) guarantees that there exists another invariant set 



£l e C W 1 for the error e, given by Jl e = {e 6 K" : |e| < a\x\, x G £l x }- 

2. Compute H and G and calculate the points (x*, e*) £ £l x x Q e and (x b , e?) £ £l x x £l e where \H\ 
and \G\ are maximized. 



3. Find a stabilizing inter-execution time r* for the linear model (4.23). Among others, one possible 
technique to compute r* was described in |Tab07] : 

r* = -* - 4f arctan ( ^ + 2 ^ \ ^ 



with: 



* = -^farctan( ai + 2 ( t° a2 ) 



B = y 4«2«o — "i 

"o = |-H"«,e*)| 

ai = \H(x* a ,e* a )\ + \G(xt,e* b )\ 

a 2 = |G(xfe,eb)| 



This time r* is a stabilizing sampling period for the original system for any initial condition lying 
in £l x (since Q x is invariant). 

4. Let r be the largest ball inside £l x , and let d be its radius. Relate the current state x with some 
point in the boundarjj^jof T via a homogeneous ray, that is, find s 6 R (the dilation in space) such 
that e rs (y) = x for some y in the boundary of T. Since we are working with the standard dilation 
and since we have an estimate r* valid for any point in the boundary of T the next execution time 
t^-(x) of the control task can be computed by using (4.19): 

T l (x)=e- p ^r* (4.25) 

4 Note that the boundary of T is an n — 1 sphere. 
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As r* can be precomputed offline, the evaluation of (4.25 ) can be performed online in a very short time. 
The estimate r^{x) represents a lower bound for the inter-execution times implicitly defined by (3.7), 



and therefore it guarantees stability for the original system. It is important to notice that the self- 



trigger technique is trying to emulate the event-trigger condition defined in (3.6). In this context, the 



conservativeness of this approach (due to the mismatch between the self-trigger and the event-triggered 



policy) relies entirely on the accuracy of the period r* guaranteeing stability for the linear system (4.23). 
That is, no conservativeness is added when the scaling law is applied. 



4.5 Example: Jet engine compressor 

To illustrate the previous results, we consider the control of a jet engine compressor. The following 
model is borrowed from [KK95 : 



-X2 



3 2 

7^ x l 



1 



T x l 



Xl = 

z z 

x 2 = 4(xi-u) (4.26) 
where xi is the mass flow, (3 is a constant positive parameter, X2 is the pressure rise and u corresponds 



to the throttle mass flow, the input for this controller. In the model (4.26) we have already translated 



the origin to the desired equilibrium point, hence the objective of the controller consists in steering both 
state variables (x\,X2) to zero. The operation region is a ball of radius 5.4 centered at the origin. A 
control law u = g(x\,x-i) is designed to render the closed loop globally asymptotically stable: 



2 

~ ^r(x\ + 1)0/ + 4v + xiV 2 ) + 2/3 2 xi 



where we have applied the nonlinear change of coordinates y 
become: 



(4.27) 

The closed loop equations 



^(x 1 2 + l)(x 1 +y) 



xi 

z 

y = -( Xl 2 + l)y (4.28) 

An ISS Lyapunov function for this system can be found using SOStools [PPSP04 : 

V = 1.46x1 - 0.35x iy + 1.16y 2 

Bounds for the derivative of V along trajectories (including measurement errors) are found using also 
SOStools: 

V < -0.74 • 10 8 |a;| 4 + 0.90 • 10 8 |x| 2 |e| 2 

where x = (xi,y) T and e = (e\,e2) T ■ Hence stability can be guaranteed if we enforce the following 
inequality: 



0.90|e| 2 < 0.74cj 2 |x| 2 



The closed loop system (4.28) is homogeneous with respect to (4.17) for £(xi,y) = -ttt- For simplicity, 

3* i ~T -1 



we pick r = 1 in the dilation vector field (4.17). We select a value of a < 1 guaranteeing stability under 
rule (3.6), for instance a = 0.33. We find an estimate of the aperiodic time sequence following the 4 



steps mentioned before. 



10 



a periodic self-trigger 



0.11 890 123 
0.22 506 68 
0.33 397 53 



Table 1: Number of executions of the control task for a simulation time of 3s. 
We define an invariant set Q, x that encloses the operation region, Q x = {x E R n |F(x) = 27.04}. 



Calculate the weighted Jacobians as defined in (4.22 ) and compute the points where the maximums 



are attained: x* = x* b = (12.5,0) and e* = e* b = (3.75,0). 

Compute the inter-execution time for the linear model: r* = 7.63ms for the selected value of a 
and the desired operation region. 



Finally, applying equation (4.251, we obtain the following formula describing the inter-execution 
times for the control task: 



T 



29 Xl (r/) +( i(x(r/)) 2 



W^»>"» ~ 5.36d(x(r/))x l( r/) 2 + d(xrf)) 2 
where d(x(r- )) is the norm of the previously measured state: 



From this formula, we can see that times tend to enlarge as the system appr oaches the equilibrium 



point, since a linear term in d(x(r- )) only appears in the denominator of (4.29). In order to show the 
effectiveness of the approach, 50 different initial conditions were considered, equally distributed along 
the boundary of the operation region. In Figures [T] and [2] we compare the behaviour of both strategies, 
periodic and self-trigger. To choose a stabilizing period for the system we select the worst case inter- 



execution time obtained from (4.29). A different way to compute a stabilizing period for nonlinear 
systems appeared in |Lai03 j ; both techniques lead to similar values for the period. The systems exhibit 
a similar behaviour under both strategies for all initial conditions tested (see Figure [T] for one particular 
initial condition, xq = (5.37,0.34)). Figure [2] shows the evolution of the input for the control system. 
At the beginning, both the periodic and self-trigger generate the same inter-execution times, but as 
the system tends to the equilibrium point the self-trigger policy increases the time between executions, 
whereas the periodic policy keeps updating the controller at the same rate. The right side of Figure [2] 
zooms the last part of the simulation, where the inter-execution times for the self-trigger strategy is 
already 24 times larger than the periodic. Hence the self-trigger implementation leads to a much smaller 
number of executions, while achieving a similar performance. The number of executions required for 
both implementations are shown in Table [T] for different values of a (and averaged over all initial 
conditions considered): the self-trigger policy reduces the number of executions by a factor of 8, for a 
simulation time of 3s. 

We also verified through simulations that the self-trigger technique is robust with respect to distur- 
bances, as expected because of the existence of an ISS Lyapunov function. Figure [3] shows the evolution 
of the states for both periodic and self-trigger strategies when sensor noise n(t) is considered (noise 
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Figure 1: Evolution of the states for self-trigger and periodic strategies. 



power being 2% of the signal power). In addition, a disturbance d(t) at the actuator level is applied at 
t = 0.7s. Again, the self-trigger strategy achieves a similar rate of decay with a much smaller number 
of executions. Figure [4] depicts the evolution of the inter-execution time under the self-trigger policy 
with and without the disturbance. When no disturbance d(t) is applied, times enlarge as the states 
approach the origin, and the inter-execution times at t = 3s rises above 210ms. For the second case, 
the disturbance steers the system far from the origin at t = 0.7s, and therefore the self-triggered task 
reduces the inter-execution times accordingly to guarantee the required performance. As the system 
reaches the equilibrium point, the inter-execution times start growing again since less executions are 
required to achieve the desired performance. Other possible real-time implementation effects that could 
degrade the performance of the system, such as jitter and input-output latency, are discussed in detail 
in [AT08 . The problem of real-time scheduling for a set of periodic, aperiodic and self-triggered tasks 
is also tackled in |AT08] . 



5 Polynomial control systems 

5.1 Scaling laws for inter-execution times of 0-related systems 

In this section we show how the results developed for homogeneous systems can also be applied to 
other classes of vector fields. In the case of homogeneous systems, the trajectories satisfy the scaling 



property (4.11). However, since the inter-execution times are determined by the evolution of the ratio 



12 




Figure 2: Control input for self-trigger and periodic strategies. 



|e|/|x|, it is sufficient to assume a scaling property just for the ratio, for instance: 

\e(t,\x , Ae )| = \e(9(x , e )t, x , e )| 
\x(t,Xx ,Xe Q )\ \x(0(xo,e o )t,xo,eo)\ 

for some function 9 : M n x ~R n — * M. and A 6 K. This observation can be exploited as follows. Let our 
open loop control system be: 

x = /(x,u) (5.2) 

A controller of the form u = k{x) is designed to render the system ISS with respect to measurement 
errors. Due to the digital implementation, the feedback control law becomes u = k{x + e) and the closed 
loop is: 

x = /(x,fc(x + e)) = /(x,e) (5.3) 
Since e = — x, the measurement error can be included in the state space representation of the system: 



x 
6 



Z(x, e) 



/(x,e) 
-/(x,e) 



(5.4) 



Let the solution to (5.4) be denoted by z, and let z = (x,e) be the state of this system. Now we 
define tt x as the projection of z of the first n coordinates (ir x o z = x) and ir e as the projection of the 
last n coordinates (ir e o z = e). We also define an output map rjz '■= | vr e 1/1^1. Then, the execution 
rule |e| = c|x| becomes rjz o z = c. We thus look for another vector field Y and an output map r/y that 
describe the evolution of times for the original system, that is ^(z(t)) = 7?y(y(t)). These ideas can be 
summarized in the following theorem. 
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Figure 3: Evolution of the states for self-trigger and periodic strategies in the presence of disturbances. 



Theorem 5.1 Consider the vector field Z : M — > TM , as defined in (5.4 ), and a vector field Y : N — > TN 
giving rise to flows z and y, respectively. If there exist maps (p : M — > N and rjy ■ N — > M suc/i that the 
following diagram commutes: 




(5.5) 



that is: 

1. TifiZ = Y o (p; 

2. rj Y o(p = rjz, 

then the inter- execution times T z (z) implicitly defined by rjz ° z(i, z) = c (with c > 0) and the inter- 
execution times T y (ip(z)) implicitly defined by rjY ° y(t, <p(z)) = c coincide: 



T -A Z ) = Ty((p(z)) 



(5.6) 



The auxiliary vector field Y is 93-related to Z, and it faithfully represents the evolution of the ratio 
Kel/I^xl of the original system. Moreover, if Y is homogeneous, scaling properties can be inferred for 
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Figure 4: Evolution of inter-execution times for the self-trigger strategy. 



the ratio |7r e |/|7r x |, and thus we can derive scaling laws for the inter-execution times of the original 
vector field Z. In the next section we describe an application of these abstract ideas to the particular 
case of polynomial control systems. 

5.2 Scaling laws for inter-execution times of polynomial systems 



As a concrete application of Theorem (5.1), we consider the class of polynomial systems. As explained 
for example in |Bai80 , any polynomial vector field can be rendered homogeneous by introducing another 
state variable we!, satisfying w = 0. For instance, the polynomial system: 

XI = X1X2 + X2 

± 2 = xi (5.7) 

can be easily rendered homogeneous of degree 1 with respect to the standard dilation by adding the 
auxiliary variable w: 

±1 = XiX 2 + X 2 W 
X 2 = XiW 

w = (5.8) 



Trajectories for the original system (5.7) can be recovered from the trajectories of the auxiliary sys- 



tem (5.8) provided that w(0) = 1. This homogenization procedure can be formalized as follows. The 
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closed loop (5.3) is assumed to be polynomial, that is: 



- ctiix-, 1 * 3 



p £ N 



>l 



r kij 6 No, 



r n ■ ■ 5i ■ ■ 



(5.9) 



ska e N °, 



VA; = 1, .., n 



1, . . . ,n 



We define I as the highest degree in any of the monomials of fi for any i = l,..,n. Each monomial 
of degree m < I is multiplied by w m ~ l . With this new state variable w we need to define a dummy 
error variable e w , since under the framework of section [4] it is necessary to have as many errors as 
states. Thus each component /, becomes a homogeneous polynomial of degree I — 1 in the variables 
xi, .., x n , w, e\, .., e n , e w . Hence the state space representation of the extended system is: 



with 



x 
w 



fi(x, w, e, e w ) — y^j— i oiijX^ 



f(x,w,e,e v 




(5.10) 



. Xn ij w l ~ mi ie 



Si ■ . 



,n 



nii 



ELi ( r kii + Sfcy ) i = maxi j mi 



(5.11) 



and w(0) = 1 so that f(x,w,e,e w ) = f(x,e). Using this procedure, the system (5.10) is homogeneous 
with respect to the standard dilation. In the context of the previous section, we can see that the vector 
field: 



f(x,w,e,e v 


-f(x,w,e,e v 




Y(x, w, e, e v 

is (/^-related to the original polynomial vector field, with ip(x, e) = (x,l, e, 0) since: 

TipZ{x, e) = 



(5.12) 



f{x,e) ' 




f(x,l,e,0) 










-f(x,e) 




-f(x,l,e,0) 











5.1 



= Y(x,l,e,0) = Yo(p(x,e) 
we define the output map r]Y(x,w,e,e w ) 



To apply Theorem 

Under this embedding, the inter-execution times for the original system (|5.3l) can be computed using 



(5.13) 

M for the extended system Y. 



Theorem 5.1, as defined in the following corollary. 



Corollary 5.2 Let t(x) be the inter- execution times for the system (5.4) implicitly defined by |e| = c|x| 
for a point x £ W 1 , and let f(x,w) be the inter- execution times for the system (5.10) under the same 
execution rule, for a point (x,w) £ M n+1 . Then, r and f are related according to: 



T (x) = X 1 - 1 f(Xx,X), VA > 



(5.14) 



with I as defined in (5.11) 
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The underlying idea consists in embedding the original system in a higher dimensional space, where 
symmetries for the flows can be established and that were not present in the original n-dimensional 



space. This corollary does not provide any scaling law for the original system (5.3), as it was the case in 
sections 4.2 and 4.4 since we always have to look at the auxiliary system (5.10) to compute r. In order 



to infer a self-trigger condition from this theorem, the same steps as in Section 4.4 can be followed, 
but now carrying the computations for a n + 1 dimensional space. Even though the extended auxiliary 
system no longer needs to be stable, existence of invariant sets (required in Section 4.4 to compute r*) 
can be established. 



Propos itio n 5.3 Consider the dynamical system (5.4), with f as defined in (5.9). IfT, is an invariant 

p2n-\ 



set for (5.4), then for any A > the set T,\ = {(Ax,A,Ae,0) G 



(x, e) G X} is an invariant set 



for the extended system (5.12) 



This proposition provides us a way to construct an invariant set where the computation of t* can 
be performed. The estimate r* of the inter-execution times can then be extended via homogeneous rays 



using Corollary 5.2 It is important to notice that, unlike for the case of homogeneous systems, it is not 
possible to guarantee that the self-trigger condition always results in longer inter-execution times than 
the periodic approach, since r* is computed in a higher dimensional space. Nonetheless, the technique 
herein explained outperforms the periodic approach for all tested examples, one of them being discussed 
in the following section. 



5.3 Example: Rigid body 

We apply the previous techniques to the control of the angular velocity for a rigid body. This example 
is borrowed from [BI89 . After a preliminary feedback and normalization with respect to the moments 
of inertia, the state space representation of such system with two inputs can be simplified to the form: 



X2 
X,3 



"1 
U2 

XiX 2 



(5.15) 



A nonlinear feedback law is designed in [BI89 to render the system globally asymptotically stable: 



U-2 



—X\X 2 - 2x 2 x 3 - X1 — X3 
2x1X2X3 + 3x§ - X2 



as proved by the Lyapunov function: 

V(x 1 ,x 2 ,x 3 ) 



(xi + x 3 ) 2 + -(x 2 -x^) 2 + x| 



(5.16) 



(5.17) 



The derivative of V along the trajectories of the closed loop system (including measurement errors) is 
bounded by: 



V < -91446|x| 4 + 147190|e| 2 |3;| 2 
where x = (xi, X2, x 3 ) T and e = {ei,e2,e 3 ) T . Hence stability will be guaranteed if: 



(5.18) 



147190H 2 < 91446cj 2 |x| 2 



(5.19) 
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for a < 1. For this particular example we select a = 0.01. The operation region is a ball around the 
origin of radius 15. We define an invariant set that encloses the operation region: 



O = {(x, e) 6 IT™ : V(x) < 25650, |e| < a\x\} 
and compute a stabilizing period for the system using the technique described in [Tab07] (equa- 



tion (4.24)). For the selected value of a and for this operation region we obtain r* = 4.5 • 10 s. 



To homogenize the original polynomial system, the auxiliary variable w is included in the state space 
representation for the closed loop: 

xi = -(xi + ei)(x 2 + e 2 )w - 2(x 2 + e 2 )(x 3 + e 3 )w - (xi + ei)w 2 
-(x 3 + e 3 )w 2 

x 2 = 2(xi + ei)(x 2 + e 2 )(x 3 + e 3 ) + 3(x 3 + e 3 ) 2 w - (x 2 + e 2 )w 2 

X 3 = X]X 2 W 

w = 



(5.20) 



This extended system is now homogeneous of degree 2 with respect to the standard dilation. In order 



to derive a self-trigger condition based on (5.14), it is necessary to find a lower bound f* for the inter- 



execution times f of the extended system. An invariant set S covering the whole operation region can 
be defined for the extended system using Proposition |5.3| 



S = {(Ax, A, Ae, 0) G M 2n+2 : [x, e) G E A , < A < 0.066} 

where S A = {(x,e) £ R 2n : V{x) < §(1 - A 2 ) + \{l - A 2 ) 2 , |e| < <j\x\}. As before, the weighted Jacobians 
for the extended system ( 5.20[ ) are computed, and the maximums are attained at: 

x* a = (1.3,2.3,-1.4), w* a = 0.066, e* = (0.031, 0.007, -0.035). 

x* h = (0.7,2.2,-1.3), 1^= 0.066, e* b = (0.028, 0.005, -0.047). 



Finally, applying the results in Corollary 5.2 the following self-trigger condition for the control law is 
found: 



T 



T 



X 



l + |x(r/)| 2 



(5.21) 



where f* = 5.1ms. Equation (5.21) shows that times tend to enlarge as the system approaches the 
origin. The rigid body system was simulated under the self-trigger and periodic implementations. 
Sensor noise was also included in the simulations (noise power being 2% of the signal power). We 
compare the behaviour of the Lyapunov functions for the periodic and self-trigger implementations in 
Figure [5] Again, both techniques achieve a similar performance. The evolution of times is displayed 
in Figure |6j At the beginning, the periodic strategy generates a larger inter-execution time than the 
self-trigger, but as the system approaches the origin, the self-trigger policy increases the inter-execution 



times according to (5.21 ). For a simulation time of 5s, the self-trigger policy reduces the total number 



of executions by a factor of 8. 



6 Appendix 



In this appendix we present the proofs of the results stated in this paper. Theorem 4.3 and Theorem 4.5 



are not proved since they are special cases of Theorem 4.11 and Theorem |4.8| respectively, for the case 
when the degree of homogeneity is constant. 
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Figure 5: Evolution of the Lyapunov function for self-trigger and periodic strategies. 

Theorem 4-8: Let X and D be vector fields on a manifold M, giving rise to flows <p :R x M —>■ M and ip : M x M 
respectively. The vector field X is homogeneous with degree function £ : M — * R with respect to D if 
and only if: 

<Pt°*P s = ip s ° 4> eP (s)t s,teR (6.1) 

with: s 

p( s ) = f £ o ^ dr (6.2) 



JO 

Proof: Let x be a point in the manifold M. Both the infinitesimal and the integral part need to 
be proved. We start with the infinitesimal part. 

4> t o ip s (x) = ip s o 4> eP(s)t (x) = ip s ° 4>( e p(s)-i) t ° M x ) 

O (f> t o V s o tp-t(x) = ip s o 4> {eP (s)_ 1)t (x) (6.3) 
We differentiate first the left hand side of the previous equation with respect to s: 

^j-<f>t o ° 4>-t{x) = T^o^^fa ^-ij) B (<j)_ t {x)) 

= T^ a0<j> _^ x )4>tDo^ s {<j)- t {x)) (6.4) 

where the last equality holds since tp satisfies ip = D{ip). Evaluating the previous expression at s = 
we obtain: 



d_ 

ds 



4> t oip s o p_ t (x) = T^_ t{x) p t D{p^ t {x)) (6.5) 

=o 
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Figure 6: Evolution of the inter-execution times for self-trigger and periodic strategies. 



Differentiating the right hand side of equation (6.3) with respect to s we obtain: 

d 



d_ 

ds 



Ms). 



-l)t 



(x) 



ds 



lps((l>(eP(.°)-l)t( x )) 
d 



(eP(s) _ 1)t (^ S ^-)-l)t(*) 

Do^ s O (f)( e p(s)_ 1)t (x) 



+ T * fB ,(.)_ 111 <*)^ X o (epW _ 1)t (x) e^k^-p(s) 



%P( s )-l)t 



ds 



Using equation (6.2), we find the derivative of p(s) with respect to s: 

d 



ds p{s) = i 



(6.6) 
(6.7) 



Since equation (6.1 ) has to be true for any value of s, it also needs to hold for s = 0, hence p(0) = 0. 

ve obtain: 

ips ° 4>( eP (s)-i)t( x ) = D ( x ) + X(x)t((x) 



Evaluating expression (6.6) at s = we obtain: 

d_ 

ds 



(6.8) 



s=0 



We now equate the right hand sides of (6.5) and (6.8), differentiate them with respect to t and finally 
evaluate them at t = 0: 



T^^tD o <j,_ t {x) = D(x) + X(x)tC(x) 



d 
dt 



t=o 



T<t,_ t ( x )(l>tD o 4>- t (x) 



d 
dt 



D(x) + X(x)t£(x) 



t=o 



[D,X}(x) = £(x)X(x) 
20 



(6.9) 



For the integral part, we need to check the other direction of the statement: 

[D,X)(x) = £(x)X(x) => <j> t o ip 8 (x) =ip s ° <j) B p{ s ) t (x) 



(6.10) 



We will integrate the first equality of Q6.10 1 with respect to t and s. First we review the definition of 
pullback of a vector field X under the mapping (see chapter 2): 



We define: 



r s X(x) = T x ^ s Xo^ s (x) 
j(s,x) = 4>* s X(x) 



and note that: 

1 (0,x) = r X(x) = X(x) 
Moreover, the function 7 satisfies the differential equation: 



d_ 
ds 



l(s,x) = ±tfX(x) = r.[D,X](x) = p a (SX)(x) 



whose solution is: 



= i>* s £,{xW s X(x) = £ o V s (x)7(s, x) 
7 (s,x) = eSo ^^ dT j(0,x) = e^ 7 (0,x) = e p ^X(x) 



according to the definition in (6.22). Hence: 

ip* s X{x) = e p ^X{x) 
The next step is to prove the following implication: 

r s X(x) = eP^X(x) => <k o ip s (x) = ip s o 4> eP{s)t {x) 



By definition of pullback, we can rewrite equation (6.16) as: 

{T x ip s )- l Xo^ s {x) = e p ^X(x) 

^Xo^ s ( x ) = T x ^ s e p ^X(x) 

Again, we define a function: 

X(t,x) = lp s o 4> eP{s)t {x) 

and show that it satisfies the following differential equation: 

j t x{t, x) = T^ p{s)tix) ij s X o </> eP(s)t (x)e^ = T^ )(( ^/ W I o </> ef>(s)t (x) 
= X o ip s o (f) eP (s )t (x) = X(x(t,x)) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 
(6.15) 
(6.16) 



(6.17) 
(6.18) 



(6.19) 



(6.20) 



So x is t ne flow of the vector field X. Since the initial conditions of x an d <p have to be the same 



(x(0,x) = ip s (x)) and the vector field X generates a unique flow, equation (6.1) is finally obtained: 

4> t o ijj s (x) = i) s o 4> eP(s)t (x) 
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Theorem J^.ll: Let x = /(x, u) be a control system for which a feedback control law u = k{x) 
has been designed, rendering the closed loop homogeneous with degree function £ with respect to the 
standard dilation vector field. The inter-execution times r : R n — > R implicitly defined by the execution 
rule |e| = c\x\, with c > 0, scale according to: 



with: 



r(e rs x) 



-P(s) 



t{x) VseM,r>0 



(6.21) 
(6.22) 



and where x € R n represents any point in the state space. 

Proof: Under the control law u = k(x), the closed loop becomes: 



x = /(x, fc(x)) 

Now we include the measurement errors as variables in the state space representation. Since e = —x 
(see (3.1)), the closed loop vector field can be written as: 

Z(x, e) = 



x 
6 



/(x,fc(x + e)) 
/(x,fe(x + e)) 



(6.23) 



The solution to (6.23) is denoted by (ftt{xo, eo) = {4>t > 01) (f° r 0o(£o, eo) = (%o, eo)), where 0f represents 
the first n coordinates of (ft and 0f represents the last n coordinates. Under this notation, the execution 
rule |e| = c|x| becomes \(ftf\ = c\(ftf\. Since / is homogeneous with respect to (4.17), the vector field (6.23) 
is homogeneous with respect to the following dilation vector field: 



D 



d ,A d 

rxi- h > rei — 

oxi oe{ 

i=l i=l 



(6.24) 



The homogeneous rays for this vector field are ift s (x,e) = (e rs x,e rs e). Since the closed loop (6.23) is 
homogeneous, cftt satisfies (4.8), i.e.: 



h(e rs x,e rs 0) =e rs eP(a)t (x,O), Vx G 



(6.25) 



Hence if we consider now the initial condition of the dynamical system to be (e rs x, e rs 0) the evolution 
of|#|M|is: 



\<ft e t(e rs x, e r 
\(ftf (e rs x, e r 



l^(.) t (^0)l 



\e rs cftl p(s)t (x,0)\ 



Therefore, the evolution of the ratio |0f|/|^f| starting at e rs (x,0) (the dilated point) is e p ^ times 
faster than when starting at (x, 0): a dilation e rs in the initial condition implies a dilation e p ^ in the 
inter-execution times according to equation (6.21). 



Theorem 5.1; Consider the vector field Z : M — > TM, as defined in (5.4), and a vector field 
Y : N — > TN giving rise to flows z and y, respectively. If there exist maps ip : M — > N and rjy '■ N — > R 
such that the following diagram commutes: 
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(6.26) 



that is: 

1. TipZ = Yo(p\ 

2. 7] Y o ip = r) Z , 

then the inter-execution times t z {z) implicitly defined by rjz ° z(t, z) = c (with c > 0) and the inter- 
execution times T y (if(z)) implicitly defined by i]y o y(t, <p(z)) = c coincide: 



Tziz) = Ty(<p(z)) 

Proof: Since Z and Y are (/9-related, the corresponding flows satisfy |Lee03] : 

ipoz(t,z) = y(t,ip(z)) 



(6.27) 



and from commutativity of Diagram 6.26 we can conclude: 



rjz o z(t, z) = i] Y o ip o z(t, z) = 7] Y o y(t, tp(z)) 



(6.28) 



Hence, since both rjz oz and rjy ° y are identical for any point z £ M, they generate the same sequence 
of times. 



Corollary 5.2: Let r(x) be the inter-execution times for the system (5.4 1 implicitly defined by 
|e| = c|x| for a point x G M n , and let f(x,w) be the inter-execution times for the system (5.12) under 



the same execution rule, for a point (x,w) G M n+1 . Then, r and f are related according to: 



r (z) = A*" 1 f(Ax,A), 



VA > 



(6.29) 



with / as defined in (5.11) 



M, then by Theorem 



5.1 



Proof: The system (5.4) is (/^-relat ed t o (5.12 ), for (p(x, e) = (x, 1, e, 0). If we define the output map 

r] Y {x,w,e,e w 



we can claim that the inter-execution times are identical: 
t(x) = f(x, 1) 



Moreover, since the auxiliary system (5.12) is homogeneous of degree / — 1, times scale according to 
Theorem 14.31 Hence: 

t[x) = f(x, 1) = A^ 1 f(Ax,A), VAX) 



Proposition 5.3 : Consider the dynamical system (5.4), with / as defined in (5.9). If X is an 
invariant set for (5.4), then for any A > the set T,\ = {(Ax, A, Ae, 0) G M 2n+2 : (x,e) G X} is an 



invariant set for the extended system (5.12). 



Proof: Let the solution of (5.4) be denoted by z(t,xo,eo) = (x, e) T . Since S is invariant, for any 



(xq, eo) G X the trajectory z stays in X, that is z(i, xq, eo) G X, for all t > 0. Now let the solution of the 
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extended system (5.12) be denoted by y(t, xq, wq, eo, e 



Ul , 



(x, w, e, e w ) T . Using the fact that y is the 



flow of a homogeneous vector field, we can conclude: 



Xz(t,x ,e ) = X 



St(t,x , l,e o ,0) 
e(t,x , l,e ,0) 



x(A' H, Xx , A, Aeo, 0) 
e(X l ~ 1 t, Xxo, A, Aeo, 0) 



(6.30) 



Moreover, we hav e w(t , Axo, A, Aeo, 0) = A, for all t > 0, as w = 0. Hence if the initial condition of the 
auxiliary system (5.12) is (Axo, A, Aeo, 0) for any (xo,eo) S S, the trajectory y stays in Sa for t > 0. 
Thus every invariant set E for the original system induces a family of invariant sets T,\ for the auxiliary 
system defined as S A = {(Ax, A, Ae, 0) G R 2n+2 : (x, e) £ £}, for any A > 0. 
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